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STRUCTURAL ANALYSIS OF WIND TURBINE ROTORS FOR 
NSF-NASA MOD-O WIND POWER SYSTEM 
by David A. Spera 
Lewis Research Center 

SUMMARY 

Preliminary estimates are presented of vibratory loads and stresses in hingeless 
and teetering rotors for the proposed NSF-NASA Mod-0 wind power system. The wind 
turbine in this system has a two-bladed rotor 38 meters (125 ft) in diameter on a 30- 
meter (100-ft) tower. System rated output is 100 kilowatts at a wind speed of 8 meters 
per second (18 mph). 

Preliminary blade design utilizes a tapered tubular aluminum spar which supports 
nonstructural aluminum ribs and skin and is joined to the rotor hub by a steel shank tube. 
Stresses in the shank of the blade are calculated for static, rated, and overload operat- 
ing conditions. Blade vibrations were limited to the fundamental flapping modes, which 
were elastic cantilever bending for hingeless rotor blades and rigid-body rotation for 
teetering rotor blades. The MOSTAB-C computer code was used to calculate aerody- 
namic and mechanical loads. 

The teetering rotor has substantial advantages over the hingeless rotor with respect 
to shank stresses, fatigue life, and tower loading. The hingeless rotor analyzed does 
not appear to be structurally stable during overloads. Therefore, until adequate reli- 
ability of associated automatic controls has been established, a teetering rotor will 
probably be required in order to achieve a long service life in a wind turbine which 
operates unattended. 


INTRODUCTION 

Large wind turbines must be designed for structural efficiency and reliability to- 
gether with minimum maintenance and weight in order to produce energy at a competi- 
tive cost. Major factors which dictate structural weight are the vibratory loads which 
act on the rotor and the tower. These unsteady loads may be aerodynamic, gravitational, 
or inertial in origin. For very large wind turbines (e. g. , diameters larger than 30 m 



(100 ft)), vibratory stresses caused by dynamic loads will probably be the governing de- 
sign consideration (ref. 1). 

This study was conducted to estimate the dynamic loads and vibratory stresses 
which may occur in two types of rotors for the proposed NSF-NASA Mod-0 wind power 
system* (ref. 2). As shown in figure 1 the wind turbine in this system will have a two- 
bladed rotor 38 meters (125 ft) in diameter on a tower 30 meters (100 ft) in height. The 
rated electrical output is to be 100 kilowatts at a wind velocity of 8 meters per second 
(18 mph). It will be located at the NASA Lewis Research Center Plum Brook Station 
near Sandusky, Ohio. 

Of particular interest to the designers of the Mod-0 wind turbine are the advantages 
of articulating (or hinging) the rotor blades, as is commonly done in helicopters. The 
simplest articulation method is "teetering, '* in which a pair of connected blades are 
fastened to the rotor shaft by means of a pivot pin. This permits a teetering (or see- 
saw) motion of the blades parallel to the shaft, which can reduce internal bending 
stresses in the blades and vibratory loads on the tower. A teetering rotor was used 
successfully in the 100-kilowatt Hutter-Allgaier wind generator system (ref. 3). The 
hub of this rotor is shown schematically in figure 2. 

In this study the advantages of a teetering rotor are determined quantitatively by 
direct comparison of vibratory loads and stresses in a teetering rotor with those in a 
hinge less design. Because of the preliminary nature of this investigation, stress anal- 
ysis was restricted to the shank area of the blades, using elementary strength- of - 
materials equations. Rotor loads on the tower were also compared for the two types of 
construction. 

A similar analysis was made by Wilcox and his coworkers together with von Karman 
on the 1250-kilowatt Smith- Putnam wind power system, in an effort to minimize vibra- 
tory loads on the tower (ref. 4). Only rigid-body motions of the rotor blades were con- 
sidered. Each blade was individually hinged and restrained by an adjustable damper. 
Their analysis indicated that the undamped condition produced much smaller vibratory 
loads than the fully damped (hinge less) condition. This conclusion was verified by their 
tests on the wind turbine. 

Two wind conditions were selected for estimating the effects of rotor teetering: the 
rated condition with a nominal wind speed of 8 meters per second (18 mph), and an over- 
load condition with a wind speed of 27 meters per second (60 mph). In both cases the 
rotor speed was 4.2 radians per second (40 rpm) with the rotor axis horizontal and 
pointing directly into the wind. A static condition was also selected as a baseline case 
for comparison purposes, with zero rotor speed and a wind speed of 27 meters per sec- 
ond (60 mph). 

* Cooperative project of the National Science Foundation and the National Aero- 
nautics and Space Administration. 
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Certain simplifying assumptions are made in this preliminary analysis. The most 
important of these is that the dominant vibration mode of the rotor blades is parallel to 
the axis of revolution. This is the so-called "flapping” mode. Vibrations in the plane 
of revolution ("lead- lag" mode) and torsional vibrations are not considered. It is as- 
sumed that lead- lag and torsional vibrations are similar for both the hingeless and 
teetering rotors and therefore can be neglected when comparing the two rotors. For the 
same reason, stress concentrations and rotor-tower dynamic coupling are also 
neglected. Because of these assumptions, loads and stresses obtained in this study are 
presented as estimates for comparison purposes only. For detailed design the effects 
of other modes of vibration, aerodynamic instabilities, rotor-tower dynamic coupling, 
and stress concentrations must be considered. 

The U.S. customary system of units was used in this analysis. Conversion to the 
International System of Units (SI) was for reporting purposes only. 


SYMBOLS 


A 

a o t0 V 

b 0 to b 3 

c 0 t0 c 4 

d 0 to 
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l 
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spar cross-sectional area 

polynomial series coefficients for spanwise distribution of weight, stiffness, 
and flapping deflection 

modulus of elasticity 

eccentricity 

force 

acceleration of gravity 
area moment of inertia 
kinetic energy function 
length of blade 
bending moment 
potential energy function 
point on spar 

blade weight per unit length 
rotor axis 
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spar radial coordinate 

spar section modulus 

station along blade axis 

tower axis 

energy 

wind velocity 

weight 

axial (flapping) displacement 
Cartesian coordinates 
flapping angle 

dimensionless weight and stiffness functions 

dimensionless spar radial coordinate 
elevation or coning angle 
pitch angle 

dimensionless axial displacement 

normal stress 

shear stress 

azimuth (constant) 

azimuth (variable) 

rotor rotational speed 

natural frequency 


PRELIMINARY ROTOR DESIGN 

Figure 3 shows the two rotor hub designs under consideration. In the hingeless 
rotor , flapping motion can occur only by elastic deformation. Bending of the blade 
changes the slope at its tip by an angle /3, which is termed the flapping angle. The 
slope at the hub is unchanged from its initial value, which is defined as the blade coning 
angle © b . In the teetering rotor the blade centerline is assumed to remain straight dur- 
ing flapping, with the entire blade rotating about the teetering hinge through the flapping 
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angle ft. At the same time the opposing blade rotates through an angle The first 
mode of flapping in this case is one of rigid-body rotation. 

The preliminary external shape of the rotor blade is shown in figure 4. The blade 
has a total length of 19 meters (62. 5 ft) and tapers from a chord of 137 centimeters 
54 in. ) at the hub to 46 centimeters (18 in. ) at the tip. The coning angle is 0. 19 radian. 
Cross sections are NACA airfoils of the 23 000 series with 0. 2 radian of twist between 
hub and tip. Blade pitch is variable, to control speed. 

The material for this preliminary blade is metal, although future blades may well 
be fabricated from nonmetallic composites. The blade structure is assumed to consist 
of a tapered aluminum tubular spar supporting nonstructural ribs and skin. The hub end 
of this spar has an outer diameter of 50. 8 centimeters (20 in. ) and an inner diameter of 
48. 3 centimeters (19 in. ). A steel tube joins the spar to the hub, forming the shank of 
the blade. Preliminary dimensions of this tube are an outer diameter of 32. 4 centime- 
ters (12. 75 in. ) and an inner diameter of 27.9 centimeters (11 in. ). 

Preliminary estimates of the weight and stiffness distributions in the metal rotor 
blade are shown in figure 5. The weight of 8900 newtons (2000 lbf) per blade is the 
upper limit for the Mod-0 wind turbine. Lighter blades are desirable because they 
would lower both dynamic and static loads, with resulting benefits in system reliability 
and cost. Substantial reductions in weight, of the order of 50 percent, may be achieved 
by using composites in future blades. No balance weights are included in this weight 
estimate. 

It is assumed that only the main spar contributes to the bending stiffness. Because 
this spar is a round tube, bending stiffnesses are equal in both the flapping and lead- 
lag directions. 


LOAD AND STRESS ANALYSIS 
Wind Turbine Coordinate Systems 

Before proceeding with the load and stress analysis it is necessary to establish 
coordinate systems for the wind turbine. These are shown in figure 6. Four coordinate 
systems are required: tower, rotor, blade, and spar. Coordinate axes in these four 
systems are designated by the subscripts t, r, b, and s, respectively. 

In the tower coordinate system the origin is at the intersection of the tower axis T 
and the rotor axis R, which establish a reference plane R-T. The x t axis is verti- 
cally upward along the tower centerline. The y t and z t axes lie in a horizontal plane, 
with z t along its intersection with the R-T plane and pointing away from the rotor. 
The wind vector is located with respect to the horizontal plane by an azimuth * w , 
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measured from the z t axis, and an elevation angle 0 . The rotor axis is located by 
the elevation angle 0 r and has an azimuth of 180°, with respect to the z t axis. 

The origin of the rotor coordinate axis is located a distance e p along the rotor 
axes from the origin of the tower axes. The axes x r and y r lie in the plane of revo- 
lution and rotate with the blades at an angular velocity Cl. The azimuth i// b locates 
Xj, with respect to the R-T plane. The z r axis is along the rotor axis, pointing to- 
ward the tower. 

The blade axis is located with respect to the plane of revolution by the azimuth angle 

and the elevation (coning) angle 0 b . The origin of the blade system is located a 
distance s along the blade axis, in a direction opposite to that of x r . The x b axis 
lies along the blade axis, pointing toward the rotor axis. The y b axis is the chord line 
of the airfoil section. It is located by the pitch angle d with respect to the plane of 
revolution and points toward the leading edge of the blade. The z b axis is usually the 
quarter-chord line. 

The spar coordinate axes are parallel to the blade axes, with their origin displaced 
a distance e g along the z b axis toward the suction side of the airfoil. Finally, any 
point p on the spar for which stresses are to be calculated is located by an azimuth 

*s« 


Blade Reactions 


The blade root reactions at the shaft centerline were obtained by using the 
MOSTAB-C computer code (ref. 5). These reactions are for one blade, either hinged or 
hingeless, vibrating only in the flapping mode under steady-state conditions. They con- 
sist of three forces and three moments referred to the rotor coordinate system: 
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These reactions must be transferred to the spar coordinate system before stresses 
can be calculated. The transfer equations are as follows: 


F x > s= F x,r cos0 b +F z,r sln0 b 

F y,s = F y,r cos 8 - (F z,r cos 0 b ' F x,r sln0 b> sin 3 
F z, s = < F Z , r cos ®b - F x, r sln ®b> cos 51 ♦ F y, r sln 8 
M x, s - M x, r cos ®b + M z, r sln ®b • e s F y, 8 
M y, s = “y, r cos a " (M z, r cos S b ' M x,r sin ®b> Bin 8 + e s F x, s - sF z, s 


M, 


= (M z> r cos 0 b - M Xj r sin 0 b ) cos 3 + M y> r sin 3 + sF^ g 


z, s 


(2a) 

(2b) 

(2c) 

(2d) 

(2e) 

(2f) 


in which 

© b coning angle 
3 pitch angle 

s distance from rotor axis to spar section 

e eccentricity of spar with respect to blade axes 
s 

The coning angle © b is assumed to be approximately equal to the average angular ele- 
vation of the blade from the plane of revolution (the x r ~y r plane). Also, aerodynamic 
and inertia loads are assumed to be negligible within the length s. 


Spar Stresses 

The spar is assumed to be a hollow round cylinder with inner radius r. and outer 
radius r . The axial normal stress <j x at a point p on the outer surface of this cyl- 
inder (fig. 6) is 


a x = 




z,s cos ^s- M y,s sin *s) 


(3a) 


in which 
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A - 4o - r l) 


spar cross-sectional area 


n f 4 a\ 

S (r Q - r. J spar section modulus 


and ij/ s is the azimuth of point p. The shearing stress r at point p is 

X, 

1 M x s r o + r o r i + r ? 

4 - -f 1 * — f 1 1 <F y>s sin * - F z> g cos *) 


The remaining components of stress are assumed to be negligible, or 

” ^r _ T x,r = 7 ip, r = ® (5) 

The axial and shearing stresses can be combined into principal and effective stresses as 
follows : 
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Rotor Loads on Tower 


The root reactions from two blades combine as follows to form the rotor loads on 
the tower, referred t° the tower coordinate system (fig. 6): 

. \ t = ( F x, 2r cos ^b + F y, 2r sin cos 0 r " F z, 2r sin 0 r (8a) 

F y, t = F y, 2r cos *b - F x, 2r sin *b (8b) 

F z, t = F z, 2r cos 0 r + (F x, 2r cos ^b + F y, 2r sin *b> sin 0 r (8c) 


M x, t = ( M x, 2r cos *b + M y, 2r sin cos °r ‘ M z, 2r sin 0 r + e r F y, t cos 0 r 


“y,t sM y,2r cos H ~ “x^r sin ^b ' e r F x,t cos 0 r " M y,t 
M z,t= M z,2r cos 0 r + e r F y,t sin0 r 


(8d) 

(8e) 

(8f) 


in which 

i{/b azimuth of first blade 

0 elevation angle of rotor axis 
r 

e r distance from origin of tower coordinates to origin of rotor coordinates 


x,2r 


= F. 


- F. 


^1 


x, r 


^ b +180° 


(9a) 


F y, 2r F y, r 


^b 


- F 


y,r 


i^ b +180 c 


(9b) 


F „ = F 
z, 2r z,r 


+ F, 


^1 


z, r 


^ b +180° 


(9c) 
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M x,2r =M x,r 


M. 


x.r 




^ b +180° 


M y, 2r M y, r 




- M 


y,r\ 


^ b +180° 


(9d) 


(9e) 


M 


z,2r 


- 




+ M 


z,r 


^ b +180° 


Of) 


and My t is the moment of the superstructure without the blades. If we assume that 
the superstructure with blades is statically balanced on the tower, 

M y, t = 2W b (e r + e b> cos 0 r (9g) 

in which e b is the axial distance along the rotor axis from the origin of the rotor coor- 
dinates to the centroid of the blades, and is the weight of one blade. 

The components of the rotor loads on the tower can be combined into a resultant 
force and a resultant moment, as follows: 


F t V F x, t + F y, t + F z, t ( 10a ) 

M t = ( iob > 


WIND LOADS 
Operating Conditions 

Loads and stresses were calculated for the teetering and hinge less rotors under 
three operating conditions: static, rated, and overload. The static condition was used 
as a reasonably severe baseline against which to compare the two dynamic conditions. 
In the static case the wind speed is 27 meters per second (60 mph) and the blades are 
stopped in the horizontal position with their chords perpendicular to the wind for maxi- 
mum aerodynamic loading. 
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In the rated operating condition the wind speed is 8 meters per second (18 mph), and 
rotor speed is 4.2 radians per second (40 rpm). Blade pitch is -0. 10 radian at the hub. 
The blades produce 150 kilowatts of mechanical power (200 hp) at this pitch. This pro- 
vides 50 kilowatts for losses and an output of 100 kilowatts of electrical power. 

In the overload condition the wind speed is 27 meters per second (60 mph), rotor 
speed remains at 4. 2 radians per second (40 rpm), and approximately 75 kilowatts 
(100 hp) are produced with -0.30- radian pitch. Under these conditions the blades are on 
the verge of a flapping instability which would cause large flapping angles and impact of 
the blade on the hub stops and possibly the tower. Flapping angle amplitude for the 
teetering blade at the overload condition is about ±10°. This value would increase 
rapidly with a slight increase in wind velocity or decrease in pitch angle. Thus, the 
overload condition represents an extremely severe operating condition which the turbine 
would experience for only a very short time, if ever. 


Wind Shear 

The earth's roughness produces a boundary layer which creates a gradient in wind 
velocity with altitude. This gradient is usually expressed in terms of the altitude to a 
fractional power. However, in the MOSTAB-C computer code the simplifying assump- 
tion is made that wind velocity V varies linearly with altitude, within the limits of the 
rotor, so that 


= V 0 (l + 6 j sin 


(10a) 


in which 

V Q nominal wind velocity at elevation of rotor axis 
6 constant which depends on ratio l/h^ 

l blade length 

hj. tower height 

s station along blade axis from rotor axis 

blade azimuth 

Defining wind shear AV as the difference between the maximum and minimum wind 
velocities, equation (10a) gives 

AV = 6V 0 


(10b) 
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It can be seen from equation (10a) that wind shear produces a periodic wind load 
with a frequency of once per revolution. This is the only frequency of blade vibration 
which can be excited by the wind shear. 

In this preliminary analysis the constant 6 was assumed to be 0. 15 for a ratio of 
blade length to tower height of 0. 62. This results in a wind shear of 1. 2 meters per 
second (2. 7 mph) at the design operating condition. This is probably not a conservative 
amount of wind shear and 5 should be increased in future computations to about 0. 20 
for high wind velocities and 0.30 for low velocities. 


Tower Wake 

An important source of periodic wind load is the aerodynamic interference created 
by the tower, which is upwind of the rotor. This is referred to as the tower wake effect. 
The magnitude of this periodic load is determined by the projected area of the tower 
structural elements, their average drag coefficient, and the sector of the rotor area 
affected by the wake. The truss tower shown in figure 1 was assumed to have a pro- 
jected area of 37 square meters (400 sq ft) when oriented diagonally to the wind. A con- 
servative value of 2. 0 was used for its average drag coefficient. Blades were assumed 
to be in the tower wake at azimuths from 345° to 15°. 

Although the tower wake has a fundamental frequency of once per revolution, it is a 
pulse load which can contain harmonics at many frequencies. Therefore, the tower 

wake is a potential excitation source for vibrations at any integer multiple of the rota- 
tional speed. 


RESULTS AND DISCUSSION 
Blade Natural Frequencies 

Vibrations of the rotor blades in service are based on their natural frequencies. In 
this analysis, only the first, or fundamental, mode of flapping vibration is considered; 
so only this natural frequency is required. Equations for calculating the flapping natural 
frequency of a rotating blade are given in the appendix. For the hingeless blade the as- 
sumed mode shape is the static deflection curve. For the teetering blade the assumed 
mode shape is a straight line. 

The variation in natural frequency with rotor speed is presented in figure 7. Both 
quantities are normalized with respect to the design speed of 4. 2 radians per second 
(40 rpm). The hingeless blade has a normalized natural frequency which varies from 
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2. 65 to 2. 88 as the rotor speed increases from zero to the design speed. Plotted in 
figure 7 are integer multiples of the normalized rotor speed, which indicate potential 
excitation frequencies. The line for three cycles per revolution (3P) crosses the natural 
frequency line for the hingeless blade at slightly under the design speed. This indicates 
that three flapping vibrations per revolution might be expected at the design speed. As 
shown in a later section, these vibrations do occur and are excited by the tower wake. 

The rigid-body motion of the teetering blade has a normalized natural frequency of 
0. 96, calculated by using equation (A10). Thus, the teetering blade experiences one 
cycle of vibration per revolution, which is excited by the wind shear, the tower wake, 
and gravity. The lead- lag frequency ratio for both the hingeless and teetering blades is 
2. 65 and does not change with rotor speed. 

These values of normalized natural frequency indicate that the metal blades are 
quite stiff and probably not prone to high vibration stresses under normal operating con- 
ditions. This can best be shown by comparing their frequencies to those of other rotor 
systems, as is done in figure 8. 

In figure 8, typical lead- lag and flapping frequency ratios are given for conventional 
propellers and several types of helicopter rotors. These data and the estimates for 
wind turbines were taken from reference 1. It can be seen that the Mod-0 blades are at 
least one integer multiple higher in normalized frequency than equivalent helicopter 
rotors or propellers. This substantially decreases the chances of exciting resonances 
or instabilities during service. 


Coning Angle 

The coning angle © b was determined to be 11° (0. 19 rad) by the procedure illus- 
trated in figure 9. This figure shows the flapping angle /3 for a single hinged blade 
without initial coning, as a function of the blade azimuth t// b - The variation is sinus- 
oidal, with a mean value of about 11° and an amplitude of 1. 5°. The mean value is de- 
termined by the direction of the resultant of the centrifugal and aerodynamic loads at the 
hinge line at the design wind condition. If this mean value is used as the coning angle 
for a two-bladed teetering rotor, root bending moments will remain essentially zero. 
One blade can experience positive flapping while the other, 180° out of phase, experi- 
ences an equal amount of negative flapping. 

The same coning angle was used for both the hingeless and teetering rotors at both 
the rated and overload wind conditions. However, 11° is not necessarily the optimum 
coning angle for operation over a range of wind conditions. Because the wind speed will 
usually be between 4 and 8 meters per second (9 and 18 mph) during operation, a com- 
promise coning angle of 7° will probably be used for the first Mod-0 blades. 
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The flapping amplitude of 1. 5° is composed of the following components: gravity 
(because of coning), 0. 8°; tower wake, 0.2°; and wind shear, 0. 5°. Each of these com 
ponents varies in an approximately linear fashion with respect to its load parameter. 
For example, doubling the wind shear constant 5 from 0. 15 to 0. 30 increases the flap- 
ping amplitude by 0. 5°. The sinusoidal variation in flapping angle is equivalent to a 
1.5 rotation of the plane of revolution about a vertical axis through the hub. Thus, the 
rotor reacts like a gyroscope to the moments caused by gravity, tower wake, and tind 
shear. These moments are all about a horizontal axis in the plane of revolution. How- 
ever, they cause a rotation about a vertical axis, analogous to gyroscope precession. 


Blade Root Reactions 

Figures 10 to 12 present the force and moment reactions at the root of one blade as 
a function of blade azimuth angle. These reactions are all referred to the rotor coordi- 
nate system and are for the rated operating condition. In these and the following figures, 

values for the teetering rotor are shown by dashed lines and those for the hingeless rotor 
by solid lines. 

As shown in figure 10 both the radial and tangential root forces for the teetering 
rotor are approximately equal to those for the hingeless rotor. The mean value of the 
tangential force F r is the aerodynamic lift on the blade in the direction of rotation. 
Superimposed on this is a sinusoidal gravity force whose amplitude is the blade weight 
The blade weight is five times the lift force, which indicates that gravity rather than 
wind is the dominant source of load at the rated condition. The mean value of the radial 
force F Xj r is the centrifugal loading on the blade root. Again, a sinusoidal gravity 
force is superimposed on it, 90° out of phase with that in the tangential direction. 

Figure 11 shows torque and twist moments at the blade root, which are also the 
same for both the teetering and hingeless rotors. The maximum weight torque occurs 
when the blade is horizontal at azimuths of 90° and 270°. This weight torque is three 
times the "working” torque required to produce the rated power of 75 kilowatts per 
blade. This emphasizes the dominant role of gravity at the design operating condition, 
e blade weight appears to be disproportionately large compared to aerodynamic forces, 
ending loads and stresses are a maximum at an azimuth of 270°, where the weight and 
working torques are in the same direction. This is the critical reaction for stress 
analysis at the rated operating condition. The twisting moment M y is primarily 

gravitational, resulting from the coning of the blade. Centrifugal untwisting is neglected 
in this analysis. 

Figure 12 shows the two root reactions which are different for the teetering and 
hingeless rotors. These are the axial force and the flapping moment, both of which are 
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affected by flapping motions of the blade. In figure 12(a) the flapping moment My r 
the teetering rotor is approximately zero because of the flapping hinge and the chosen 
coning angle. However, the hingeless blade experiences a significant oscillatory flap- 
ping moment with three peaks per revolution. This is the 3P elastic vibration indicated 
by the natural frequency chart (fig. 7). 

The maximum flapping moment in the hingeless blade occurs at an azimuth of 30 . 
This indicates that maximum deflections will lag maximum loads by a phase angle of 
about 30° in this blade. The maximum flapping moment is composed of gravity, wind 
shear, and tower wake components in the following proportions: gravity (IP), 25 per- 
cent; wind shear (IP), 20 percent; and tower wake (3P), 55 percent. The tower wake 
pulse apparently contains a significantly large third harmonic, which excites the hinge- 
less blade at its natural frequency. 

According to equation (9e) the resultant load on the tower from two blades is the 
difference of flapping moments 180° out of phase in azimuth. As shown by points 1 and 
2 in figure 12(a) the maximum and minimum flapping moments are also 180° out of phase 
in the hingeless rotor. Thus, the tower will be subjected to a peak bending moment 
which is the sum of the two most severe flapping moments for one blade. 

In figure 12(b) the axial force on the hingeless rotor blade also shows a 3P vibration 
about the average aerodynamic drag force. On the other hand the teetering blade ex- 
hibits a smooth sinusoidal variation in axial force at its root, which is caused by inertia 
forces accompanying its rigid-body vibration. 

In summary, four of the six reactions at the blade root are the same for both the 
teetering and hingeless rotors at the design operating condition. Included in these is the 
critical reaction for stress analysis, which is the torque moment. Gravity loads are 
dominant over aerodynamic and inertia loads. However, the hingeless blade experiences 
a 3P elastic flapping vibration, which could result in tower excitation or high vibratory 
blade stresses at higher wind speeds. 


Blade Shank Stresses 

The axial stress ct x is the dominant stress component in the effective stress 
(eq. (7)). The shear stress r x ^ is small in the blade shank area. Figure 13 shows 
the critical axial stress cycles at two locations in the shank of the blade under the 
design operating conditions. These stress cycles would be the determining factor in the 
fatigue life of the blade root and probably of the entire blade. The absolute values of 
stress shown here are not significant except for the fact that they are low enough to in- 
dicate that the preliminary blade design is feasible for both types of rotors. These 
stresses do not include a stress concentration factor. 
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Point A is 1.0 meter from the rotor axis on the trailing edge of the steel tube which 
joins the airfoil section of the blade to the hub. Preliminary dimensions of this tube are 
an outer diameter of 32. 4 centimeters (12. 75 in. ) and an inner diameter of 27. 9 centi- 
meters (11. 0 in. ). 

Point B is 1.5 meters from the rotor axis on the trailing edge of the tapered alumi- 
num spar where it joins the steel shank tube. Preliminary dimensions of the spar at 
this radial station are an outer diameter of 50. 8 centimeters (20. 0 in. ) and an inner 
diameter of 48. 3 centimeters (19. 0 in. ). 

As shown in figure 13, critical stresses are approximately equal for the teetering 
and hingeless rotors. The reason for this equality is that the critical root reaction at 
the design condition is the torque moment shown in figure 10. This moment is the same 
for both rotors because it is determined by the rated power and the blade weight. 

Critical stress cycles for the static, rated, and overload conditions are summarized 
in table I. These stresses are expressed in terms of baseline stresses a A and 
obtained for the static condition. In this condition the blades are stationary in the hori- 
zontal position with their chord planes facing the wind. The weight of the blade and the 
aerodynamic drag produce the baseline stresses in the table. 

As was shown in figure 13 the critical stress cycles for the teetering and hingeless 
rotors are approximately equal at the rated condition. However, at the severe overload 
condition the stresses in the teetering blade are only about one-half those in the hinge- 
less blade. Moreover, the cyclic frequency for the teetering blade is only one-third that 
for the hinge less blade. 

During overload conditions the stress and frequency advantages of the teetering 
blade over the hingeless blade combine to greatly reduce its relative rate of fatigue 
damage. As a preliminary estimate the damage rate in the teetering blade during over- 
load would probably be less than one-twentieth that in the hinge less blade. In other 
words, if the wind turbine were to approach an instability condition in which the blades 
begin moderate to large flapping vibrations, fatigue damage to a teetering blade would 
be minor compared to that in a hingeless blade for equal periods of time. 

An estimate of the relative importance of fatigue damage at the rated condition to 
that at overload can also be made from the data in table I. The overload stresses are 
considerably larger than the rated stresses for both types of rotors. This means that 
the blades must be designed for infinite life at the rated condition in order to have a 
reasonable cumulative fatigue life. Thus, the overall fatigue life of the rotor would not 
be determined by the blade stresses at the rated condition but by the stresses during 
gusts, control malfunctions, very high winds, and other abnormal operating conditions. 
On this basis, the life of the teetering rotor would be substantially longer than that of 
the hingeless rotor. A spectrum of loads from rated to overload is required for a quan- 
titative estimate of the fatigue life advantage of the teetering rotor. Qualitatively, it is 
probably a factor of at least 5 to 1. 
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Rotor Loads on Tower 


In figures 14 and 15 the loads exerted on the tower by the teetering and hingeless 
rotors are compared for the design condition. These loads are referred to the tower 
coordinate system, which is nonrotating (fig. 8). As shown in figure 14 the average 
forces on the tower are similar for the two rotors. However, the vibratory components 
of these forces are somewhat larger for the teetering rotor. This is particularly true 
of the side force Fy t . 

Moment loads are compared in figure 15. The shaft torque M z t is identical for 
both rotors. As shown in the figure this torque experiences a drop in magnitude of about 
20 percent each time a blade passes through the tower wake at azimuth angles of 0° and 
180°. Twisting moments on the tower M x ^ are somewhat larger in amplitude and 
twice as frequent for the hingeless rotor. The pitching moments on the tower My ^ are 
considerably different for the two types of rotors. At an azimuth of about 30° the hinge- 
less rotor will exert a large pitching moment. This moment results from the two blade 
reactions shown in figure 12 as points 1 and 2. This moment peak is the critical rotor 
load on the tower. 

For comparison purposes, the maximum values of the resultant force load and the 
resultant moment load were calculated by means of equations (10). These are summar- 
ized in table II for the three operating conditions: static, rated, and overload. 

The resultant static force F + _ is a maximum when the stationary blades are hori- 

s 

zontal and the chord planes are perpendicular to the wind. The maximum resultant 
static moment Mj. occurs with the blades vertical. The lower blade is assumed to be 
completely shielded from the wind by the tower. 

At the rated wind condition the maximum force load is the same for both rotors and 
is about the same as the static load at the higher wind speed. Vibratory components are 
small and are neglected. Rated moment loads are one-third to one-half the static values 
for teetering and hingeless rotors, respectively. 

The overload condition produces a substantial increase in tower loads for both 
rotors. Most significant is the maximum resultant moment load of 4. 80 M t g from 
the hinge less rotor. Under the same condition the teetering rotor exerts only one- 
fourth the moment on the tower. 

Thus, as for the blade stresses, there is little difference between the two rotors as 
measured by their loads on the tower at the rated condition. However, at the overload 
condition, loads from a teetering rotor will be much less severe than those from a 
hingeless rotor. 
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CONCLUSIONS 


Estimates were made of the vibratory loads and stresses in two types of rotors for 
the NSF-NASA Mod-0 wind turbine. Both hingeless and teetering rotors were analyzed 
under static, rated, and overload conditions. Blade vibrations were limited to the flap- 
ping mode. The following major results and conclusions were obtained: 

1. A teetering rotor will probably be required in order to achieve a long service 
life in a wind turbine which operates unattended and is exposed to periodic overload con- 
ditions. This type of rotor has substantial advantages over a hingeless rotor with re- 
spect to shank stresses, fatigue life, and tower loading. 

2. The hinge less rotor analyzed does not appear to be structurally stable during 
overloads which could result from high winds, severe gusts, or control malfunctions. 
Therefore, unattended operation of a large wind turbine with a hingeless rotor is not 
recommended until high reliability of associated automatic controls has been established. 

3. If operations do not exceed rated conditions in severity, stresses and tower loads 
are about the same for both the teetering and hingeless motors analyzed in this study. 

4. For a spectrum of operation between rated and overload the teetering rotor will 
probably have at least five times more fatigue life than the hinge less rotor. 

5. For equal structural stability the tower for a hingeless rotor would have to be 
heavier and stiffer than that for a teetering rotor, 

6. The Mod-0 rotor, either hinge less or teetering, will probably be substantially 
less susceptible to dynamic instability than helicopter rotors or propellers. 

7. The main source of loads and stresses in both rotors is the blade weight, which 
is disproportionately large compared to aerodynamic forces. Efforts should be made to 
substantially reduce the weight of the blades, possibly by using nonmetallic composite 
materials. 

Loads and stresses obtained in this study are presented as estimates for compari- 
son purposes only. For detailed design the effects of other modes of vibration, aerody- 
namic instabilities, rotor-tower dynamic coupling, and stress concentrations must be 
considered. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, December 4, 1974, 

778-00. 
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APPENDIX - EQUATIONS FOR NATURAL FREQUENCY 


OF BLADES IN FLAPPING MODE 

The derivation of the following equations follows the general procedures of refer- 
ences 6 and 7, with the effect of coning angle added. If we assume that the vibrating 
blade does not dissipate energy, its maximum potential energy equals its maximum 
kinetic energy, or 


U = U, 

p, max k, max 


(Al) 


in which p and k denote potential and kinetic energies, respectively. The maximum 
potential energy of flapping vibrations is as follows, for a blade with its root at the shaft 
centerline: 


U, 


p, max 




dx 




qw^ dx + W t w^ 


(A2) 


All coordinates are in the blade system, and 

l length of blade 

E modulus of elasticity 

I moment of inertia of a cross section about its y (lead- lag) axis 

yy 

x axial coordinate 

w maximum flapping displacement 

O rotational speed 

0^ coning angle 
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q weight per unit length 
W|. concentrated weight at tip 

w t maximum tip displacement 

If we assume the motion to be simple harmonic, the maximum kinetic energy of the 
blade is 


U, 


k, max 


- 4 ( 7 ' 

2 g Fo 


qw 2 dx + W t w 2 


(A3) 


in which o> w is the unknown natural flapping frequency. The following dimensionless 
ratios are now defined: 


V = x/l 

(A4a) 

11 

* 

\ 

(A4b) 

n - E V w b' 2 

(A4c) 

y 2 = q lpN h 

(A4d) 

r x =fi 2 Z/g 

(A4e) 

r 2 = w t /w b 

(A4f) 


in which W b is the total weight of one blade. Introducing equations (A4) into (A2) and 
(A3) then gives the potential energy equation 


W b* 

U P 

p, max 2 


(A5a) 


in which 




Vl<«") 2 + rJy 2 

V cos 2 0 b J 

V o 9 9 

3 UT dr? - r sin^0 b 

+ r 2 (^>) 2 cos 2 0 b l 



- 

J_ 


d v 


r i r 2^t sin20 b (A5b) 
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The kinetic energy equation becomes 



max 



w/ 

2g 


K 


(A6a) 


in which 


K = ^ y 2 i 2 drj + T (A6b) 

(Primes denote differentiation with respect to rj . ) Combining equations (Al), (A4a), and 
(A5a) yields the natural flapping frequency of the rotor blade as 


CO 


W 



(A7) 


Hinged Rotor Blade 

The first flapping mode of a hinged blade is one of rigid-body rotation. This is also 
the case for a teetering rotor, providing that the proper coning angle is selected to make 
the average flapping angle equal to zero. The mode shape and its derivatives then be- 
come 


f; = 

V =0 > 

V I = o. 


(A8) 


in which /3 is the flapping angle. 

Combining equations (A5b), (A6b), and (A8) gives 


P=|3 Tj cos 2© b M 72 1 ? + T g 




K 




y 2 n dTj + r 2 


(A9a) 

(A9b) 
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Substituting equations (A 9) into (A7) then gives the natural flapping frequency for the 
hinged blade: 


u>. 


w 



cos 2© b 


or 


<^ w = cos 20 b 

An equation similar to (A10) is given in reference 4. 


(A10) 


Hingeless Rotor Blade 


The first flapping mode shape of the hingeless blade can be approximated by its 
static deflection curve. By Rayleigh’s principle, this approximation will lead to a cal- 
culated natural frequency which is slightly higher than the actual frequency. Polynomial 
series can be used to express the static deflection curve and its derivatives in terms of 
the distribution of stiffness and weight along the blade as follows: 


4 = a Q + ajTj + a 2 7} 2 + a 3 ?7 3 + a 4 ^ 4 + a^ 5 

O O A 

4’ = a^ + 2a 2 Tj + 3a 3 r/ + 4a + 5a > 
4” = 2a 2 + 6a 3 ?] + 12a^?j 2 + 20a.gq^ 


(All) 


For a hingeless blade with zero deflection and rotation at the shaft centerline, the poly- 
nomial coefficients are as follows : 
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in which b Q to bg 
as follows: 


a 0 ~ a l " o 


a 2 " 


2br 


a 3 - 7T <°1 - 2a 2 b l> 
6b 0 


— ^2 " ^ 2^2 “ ^9.gb|) 

12b Q 


a 5 - — — ( c 3 " ^agbg - 6agbg - 12a ^bj) 


r 


20b 


j 


are curve-fit constants for the dimensionless stiffness 


and 


El, 


yy _ 


w b r 


= ?1 * b 0 + M + b 2^ + V*' 



c 0 - "( c i + c 2 + c 3 + c 4^ J 


(A 12) 


distribution 


(A 13) 


(A 14) 


The constants <3 q, dp and dg are curve- fit constants for the dimensionless weight 
distribution, as follows: 
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- (A15) 


-SLL = 

Wu 


y 2 = d 0 + d l ?} + d 2^ 


The natural frequency of the hingeless blade is then obtained by substituting equa- 
tions (All) to (A15) into (A5b) and (A6b). 
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TABLE I. - SUMMARY OF BLADE ROOT STRESSES 


[stress concentration factor, 1.00.] 


Radial 

station 

Operating 

condition 

Wind speed 

Axial stress cycle and frequency a 

m 

ft 

m/sec 

mph 


Hinge less rotor 

b 1.0 

3.3 

Static 

27 

60 

m 

1.00 <J A 

Rated 

■ 

18 


(0. 20±0. 42) <J A ; 
0. 67 Hz 

Overload 


60 


(0.16*2.24) a A ; 
2.0 Hz 

c l. 5 

5.0 

Static 


60 


1. 00 CTg 

Rated 

■ 

18 

(0.25±0. 41) u n ; 
0. 67 Hz 

(0. 25*0. 41) ct b ; 
0. 67 Hz 

Overload 

Q 

60 

(0. 15±1. 32) ff B ; 
0.67 Hz 

(0.25*2.32) or B ; 
2.0 Hz 


a c A = 8600 N/cm 2 (12 500 lbf/in. 2 ); <r 0 = 4900 N/cm 2 (7100 lbf/in. 2 ). 
"Point A, where steel shank tube joins hub. 
c Point B, where tapered aluminum spar joins steel shank tube. 


TABLE n. - SUMMARY OF ROTOR LOADS ON TOWER 


Type of 
load 

Operating 

condition 

Wind speed 

Maximum resultant load and critical component frequency 51 

m/sec 

mph 

Teetering rotor 

Hinge less rotor 

Force 

Static 

27 

60 

1 - 00F t,s 

l -°° \s 

Rated 

8 

18 

°-" F t,s 

°- 99F t,s 

Overload 

27 

60 

1.48 F t g ; 1.34 Hz 

1.59 F t 2.67 Hz 

1 , s 

Moment 
^ 

Static 

27 

60 

1.00 M^ 

1. 00 M^ g 

Rated 

8 

18 

0.35 Mj. g ; 1.34 Hz 

0.48 L34 Hz 

Overload 

r» n ortn it /m 

'27 

60 

1.15 M t g ; 1.34 Hz 

4.80 M t 2.67 Hz 


F t, s = 38 200 N (860 ° lbf ); M t, s = 132 000 N ‘ m < 97 500 















Figure 1. • NSF-NASA MOD-O wind power system. Rated power output, 100 kilowatts; rated wind speed, 8 meters per second 
(18 mph). 
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Figure 3. - Schematic illustrations of hingeless and teetering wind turbine rotors. 
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Figure 8. - Blade frequencies for various rotor systems. (Data from ref. 1. 
except NS F- NASA MOD- a I 



Figure 9. - Flapping motion of a hinged Mod-0 blade at rated wind speed of 
8 meters per second (18 mph). 
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Figure 12. - Axial force and flapping moment at a blade root at rated wind speed of 
8 meters per second (18 mph). 
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Blade azimuth, deg 

Figure 15, - Moments on tower from rotor at rated wind speed of 8 meters per second (18 mph). 
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